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ABSTRACT 

The planet orbiting r Boo at a separation of 0.046 AU could produce a reflected light flux as bright as 
1 x 10~ 4 relative to that of the star. A spectrum of the system will contain a reflected light component 
which varies in amplitude and Doppler-shift as the planet orbits the star. Assuming the secondary 
spectrum is primarily the reflected stellar spectrum, we can limit the relative reflected light flux to be 
less than 5 x 10 -5 . This implies an upper limit of 0.3 for the planetary geometric albedo near 480 nm, 
assuming a planetary radius of 1.2 i?j up - This albedo is significantly less than that of any of the giant 
planets of the solar system, and is not consistent with certain published theoretical predictions. 

Subject headings: planetary systems — stars: individual [t Bootis) — techniques: spectroscopic 



1. INTRODUCTION 

Radial velocity surveys of nearby F, G, K and M dwarf 
stars have revealed eight planets (Mayor & Queloz 1995; 
Butler et al. 1997; Butler et al. 1998; Fischer et al. 1999; 
Mayor et al. 1999; Queloz et al. 1999) which orbit their 
parent stars with a separation of a <J 0.1 AU. These close- 
in extrasolar giant planets (CEGPs) may be directly de- 
tectable by their reflected light, due to the proximity of the 
planet to the illuminating star. In this Letter, we present 
the results of a spectroscopic search for the reflected light 
component from the planet orbiting the star r Boo. The 
motivation to attempt such a detection for a CEGP is 
strong: It would constitute the first direct detection of a 
planet orbiting another star. It would yield the orbital 
inclination, and hence the planetary mass, and would also 
measure a combination of the planetary radius and albedo, 
from which a minimum radius can be deduced. Further- 
more, it would open the way to direct investigation of the 
spectrum of the planet itself. Conversely, a low enough 
upper limit would provide useful constraints on the radius 
and albedo of the CEGP. 



2. REFLECTED LIGHT 

2.1. Photometric Variations 

In order to calculate the predicted flux ratio of the 
planet relative to the star, let R p denote the planetary 
radius, R s the stellar radius, a the physical separation, 
and a the angle between the star and the observer as seen 
from the planet. The observationally useful quantity is the 
geometric albedo p, which is the flux from the planet at 
a = divided by the flux that would be measured from 
a Lambert law (ie. perfectly diffusing; see, for example, 
Sobolev 1975) disk of the same diameter, located at the 
distance of the planet. In the case that R p <C R s <C a, the 
ratio e of the observed flux from the planet at a = to 



that of the star is 
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The value of p depends on the amplitude and angular 
dependence of the various sources of scattering in the 
planetary atmosphere, integrated over the surface of the 
sphere. For a Lambert law sphere, p = 2/3, whereas 
for a semi-infinite purely Rayleigh scattering atmosphere, 
p = 3/4. The geometric albedos at 480 nm of Jupiter, Sat- 
urn, Uranus and Neptune are 0.46, 0.39, 0.60, and 0.58, 
respectively (Karkoschka 1994). 

We treat the orbit as circular, since the observed orbit 
of t Boo has an eccentricity less than 0.02 (Butler et al. 
1997). We neglect occultations, since a transit would pro- 
duce a ~ 0.01 mag photometric dip, and is ruled out by 
Baliunas et al. (1997). We take the orbital phase $ G [0, 1] 
to be at the time of maximum radial velocity of the star. 
The phase angle a £ [0, tt] is then defined by 

cos a = — sin i sin 2ir$ (2) 

where i € [0,7r/2] is the orbital inclination. The flux from 
the planet at a phase angle a relative to that at a = is 
denoted by the phase function (f>(a). In the case of a Lam- 
bert law sphere, the phase-dependent flux ratio f(Q, i) is 
given by (Sobolev 1975) 



= e<j)(a) =p 
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sin a + (tt — a) cos a 



(3) 

For this analysis, we assume the phase variation of the 
reflected light is described by equation 0. The phase func- 
tions of the gas giants of our solar system are well approx- 
imated as Lambert spheres (see, for example, Pollack et 
al. 1986). 

In the case of r Boo, Baliunas et al. (1997) can exclude 
a sinusoidal photometric variation at the planetary orbital 



2 



period with a peak-to-peak amplitude of 0.4 millimag or 
greater. The predicted variation due to a highly reflective 
companion of Jupiter size is ~ 0.1 millimag. If proposed 
photometric satellite missions (Matthews 1997; Rouan et 
al. 1997) can achieve a precision of ~ 10 /xmag with sta- 
bility over timescales of a few days, they could measure 
this photometric modulation, as discussed by Charbon- 
neau (1999a). 

2.2. Spectroscopic Variations 

We assume that r Boo has a stellar mass of M s = 
1.2Mq, based on its spectral classification as an F7 V 
star (Perrin et al. 1977). It has a [B — V) color of 
0.48, consistent with the spectral classification. The ra- 
dial velocity observations (Butler et al. 1997; Marcy 
1997) provide the orbital period (P = 3.3125 d), phase 
(T$ =0 = 2450526.916 JD), eccentricity (e = 0) and am- 
plitude (K s — 468 m s _1 ), from which a semi-major 
axis of a = 0.0462 AU and a planetary mass of M p = 
3.89-Mj U p/ shn are calculated. The radial velocity of the 
planet relative to that of the star is 

u p ($, i) = -K s Mfl + M P C os 2tt$. (4) 

IVlp 

This has a maximum amplitude of |fp(3>, i)\ — 152 kms -1 . 

Thus the spectrum of the system could contain a sec- 
ondary component which varies in amplitude according to 
equation || and in Doppler-shift according to equation || 
Charbonneau, Jha & Noyes (1998) demonstrate that the 
effect of the reflected light component on the line profile 
bisector is not far from current observational limits for a 
CEGP of high reflectivity. This is an alternate technique 
which may be used to directly detect or limit the reflected 
light from a CEGP. 

2.3. Tidal Locking Effects 

Baliunas et al. (1997) use the activity-rotation relation 
of Noyes et al. (1984) and the mean Ca II flux of r Boo 
to predict a stellar rotation period of 5.1 days. Analysis 
of the observations of the Ca II H & K lines by Baliunas 
et. al (1997) yield a weakly detected period of 3.3 ± 0.5 d, 
consistent with the observed orbital period of the planet, 
implying that the star and planet form a tidally locked 
system. Marcy et al. (1997) demonstrate that, in the case 
of r Boo, a convective envelope of mass Mqe ~ O.OIMq 
could become tidally locked in less than the age of the 
system. If so, then there is no relative motion of any 
point on the surface of the planet relative to any point 
on the surface of the star. In this case, the planet re- 
flects a non-rotationally-broadened stellar spectrum, with 
a typical line width dominated by the stellar photospheric 
convective motions (~ 4 kms" 1 ; Baliunas et al. 1997). 
Thus, one might expect relatively narrow planetary lines 
superimposed on much broader stellar lines. 

3. TARGET SELECTION AND OBSERVATIONS 

Several considerations entered into the choice of r Boo 
(HR 5185, HD 120136) as the optimal candidate for this 
experiment. Firstly, the semi-major axis of r Boo was 
smaller than that of the other three CEGPs (51 Peg, v 
And, & p 1 Cnc) known at the time, which is desirable since 



the relative amplitude of the reflected light decreases with 
the square of the planet-star distance. Secondly, the vi- 
sual brightness of r Boo is greater than either 51 Peg or p 1 
Cnc. The photon noise of the star is the dominant source 
of noise in the experiment, and a brighter star allows for 
a more precise determination of the stellar flux in a given 
amount of observing time. Thirdly, as discussed above, 
it may be that the star is rotating with the planetary or- 
bital period. If so, the planetary spectral features would 
be much sharper and deeper than those of the primary, 
which might facilitate their separation. 

We observed r Boo for three nights (1997 March 20 to 
March 22) using the HIRES cchelle spectrograph (Vogt et 
al. 1994) mounted on the Keck-1 Ten-Meter Telescope at 
the W. M. Keck Observatory located atop Mauna Kea in 
Hawaii. These nights were carefully selected based on the 
phase of the companion's orbit. The spectral range used 
in this analysis was 465.8 nm to 498.7 nm, and the ob- 
servations were made at a resolution R = X/5X of either 
60 000 (March 20) or 45 000 (March 21 & 22). 

Since the apparent magnitude of r Boo is 4.5 mag, the 
high flux from the star would saturate the detector pix- 
els for an exposure time less than the readout time of 
the CCD. To avoid this readout-time-limited scenario, the 
cross-disperser was slowly trailed during each observation 
so as to spread the photons over roughly 30 pixels. This 
allowed for typical exposure times of 300 seconds, which 
resulted in a duty cycle of ~ 70%. In all, 154 spectra of 
r Boo were obtained, with a nightly breakdown of 38 for 
March 20, 32 for March 21 and 84 for March 22. Cloudy 
weather degraded the number and quality of the spectra 
on March 21. 

4. DATA ANALYSIS 
4.1. Extraction of the One- Dimensional Spectra 

Since the extraction of the one-dimensional spectra 
from the two-dimensional exposures must be accomplished 
without introducing systematic errors above the level of 
1 x 10~ 4 per dispersion element, it was necessary to create 
an entirely new and independent set of extraction codes 
specific to this experiment. By so doing, we were able to 
treat the sources of systematic noise particular to the Keck 
HIRES and these observations, as well as have the neces- 
sary control in identifying sources of contamination as our 
understanding of the data proceeded. 

To extract the one-dimensional spectra from an indi- 
vidual frame, the following algorithm was applied: The 
bias was subtracted and the non-linear gain was cor- 
rected. A two-dimensional scattered light model was de- 
rived by fitting the inter-order scattered light, and sub- 
tracted. The two-dimensional flat-field correction was ap- 
plied, and the orders were extracted by summing along the 
cross-dispersion direction, making use of windows which 
we had produced to identify the location of both the spec- 
tral orders and the cosmetic defects from internal reflec- 
tions and a felt-tip pen mark. The one-dimensional spec- 
tra were then corrected for cosmic rays by cubic spline in- 
terpolation across contaminated pixels. A low amplitude 
source of high frequency noise in the extracted spectra was 
corrected for by applying a narrow notch filter in Fourier 
space. The typical signal-to-noise-ratio (SNR) was ~ 1500 
per dispersion element. The wavelength solution was de- 
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rived from extracted Th & Ar emission line spectra taken 
throughout the observing run. 

4.2. The Model 

The model is that the data contain a secondary spec- 
trum, identical to that of the primary, but Doppler-shifted 
due to the orbital motion of the planet and varying in am- 
plitude with the angle subtended between the star, the 
planet and the observer. The key to the method is to first 
produce a stellar template spectrum, and then make use 
of the orbital parameters from the radial velocity obser- 
vations to calculate a model for a given observation taken 
at a particular orbital phase. The methods we briefly de- 
scribe here will be presented in detail in an upcoming pa- 
per (Charbonneau 1999b). 

The high SNR stellar template spectrum was produced 
for each of the two instrumental resolutions by combining 
all of the extracted spectra. Initially, a high SNR spectrum 
was chosen and an optimized model was found which cor- 
rected each observation to this reference (allowing for vari- 
ations in the wavelength solution and instrumental profile, 
and low- frequency spatial variations of the continuum) . A 
summed stellar template spectrum was produced, and this 
process was iterated twice, beyond which point the errors 
were no longer significantly reduced by further iteration. 
The errors in the summed stellar template were 1.2 times 
the expectation from photon noise, indicating a precision 
of ~ 1 x 10~ 4 per dispersion element. We note that this 
may well comprise the most precise visible stellar spectrum 
for any star other than the Sun. 

For each observed spectrum, we first modify the stellar 
template spectrum in order to correct for the aforemen- 
tioned variations in the wavelength solution and instru- 
mental profile, and low-frequency spatial variations of the 
continuum. Note that we wish to interpolate the stellar 
template spectrum, and not perform the reverse proce- 
dure and interpolate the observed spectra, since the stellar 
template spectrum is at a much higher SNR. Then, if we 
denote by S this modified stellar template spectrum and 
by A the wavelength solution, the model at a given pixel 
j of an observed spectrum taken at an orbital phase $ is 
described by 



r Boo and includes the instrumental effects. Thus we also 
investigated the model 



M j = 



S(Xj) + e </>($,*) s(Xj 
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(5) 



The two unknown parameters are {e,i}- The situation in 
which there in no reflected light from the planet is equiv- 
alent to e = 0. In this case, the observed spectra are best 
fit as replicated stellar spectra. 

As noted earlier, the stellar rotation period may be 
the orbital period of the planet, and hence the reflected 
spectrum may be composed of non-rotationally-broadened 
lines. The instrumental resolution will smear all spec- 
tral lines to a width of ^ 7 kms -1 . Several exposures of 
the sharp-lined F8 V star 36 UMa (HR 4112, HD 90839, 
V = 4.84, B — V — 0.52) were combined and corrected to 
the Doppler-shift of r Boo to produce a stellar template 
spectrum, S'. The spectral differences between an F7 V 
and an F8 V star are insignificant for the purposes of this 
analysis. The spectrum of 36 UMa serves as an excellent 
mock-up for the non-rotationally-broadened spectrum of 
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1 + e <£($,i) 



(6) 



where 7 is a normalization factor. 

The model was evaluated by calculating the x 2 parame- 
ter as a function of {e, i}. The minimum x m in i s subtracted 
off to define Ax 2 = X 2 — Xmin- The confidence levels in the 
allowed values of the parameters are described by draw- 
ing contours of fixed A% 2 at a desired set of significance 
levels. The confidence levels were tested for a given choice 
of {e,i} by directly injecting a reflected light secondary 
at the correct amplitude and Doppler-shift into each ob- 
served spectrum. At e > 10 -3 and high inclination, the 
secondary can be detected at the 99% confidence level with 
only one spectrum. At e ~ 10~ 4 , the planet is recovered 
only by considering all of the spectra, and the uncertainty 
in the parameters is significantly greater. Tests showed 
that the planet could be recovered for i ;> 10°. 

A second test was provided by the detection of solar 
contamination employing a model similar to the one de- 
scribed in equation H, but with the modification that the 
secondary spectrum is at a constant (but unknown) ampli- 
tude and Doppler-shift. Solar contamination was detected 
at the Doppler-shift between the Sun and r Boo, and at a 
relative amplitude of 10~ 3 , in the spectra taken on March 
21. The source of this contamination was reflection of 
the solar spectrum off the Moon and subsequently off the 
clouds which were present throughout the night. The ex- 
clusion of the contaminated spectra from the reflected light 
analysis did not greatly reduce the statistical significance 
as these spectra contained only 10% of the photons of the 
entire data set. 

5. RESULTS AND DISCUSSION 

We find no evidence for a highly reflective planet orbit- 
ing r Boo. For i ;> 10°, we can constrain the reflected 
flux ratio e <; 8 x 10 -5 at the 99% confidence level, un- 
der the assumptions that the reflected light spectrum is a 
copy of the stellar spectrum. For i ;> 70°, this improves to 
e < 5 x 10~ 5 . Assuming a planetary radius of 1.2 i?j up 
(Guillot et al. 1996), this limits the geometric albedo 
to p < 0.3. Figure 1 shows the precise limit of the re- 
flected light amplitude as a function of orbital inclination. 
Under the assumption that the secondary reflects a non- 
rotationally-broadened version of the stellar spectrum, this 
limit becomes stronger for high inclinations. The particu- 
lar shape of a given confidence level in Figure 1 results from 
the interplay of the orbital phases and statistical weights 
of the set of spectra. The upper limit imposed is set pri- 
marily by the last night of observing (March 22), when 
the planet was near a phase of $ = 0. The dip down 
to stronger constraints on the flux ratio at an inclination 
i ~ 15° results from the first night of observing (March 
20) when the planet was near inferior conjunction: Only if 
the planet is at low inclinations will it be expected to con- 
tribute a reasonable reflected light signal and hence allow 
us to significantly differentiate between models. 

At very low inclinations (i < 10°), this experiment is 
not able to exclude even very bright companions due to 
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both the lack of a significant Doppler-shift between the 
primary and the secondary, and the lack of a phase vari- 
ation in the light from the secondary. However, these low 
inclination orbits may be excluded under a further consid- 
eration: If the axis defined by the stellar rotation is the 
same as that of the orbit of the planet, then the observed 
v sin i ~ 15 kms -1 for the star would imply a true rota- 
tional velocity of greater than 50 kms -1 for i <^ 17°. Such 
high rotational velocities are not observed (Gray 1982) for 
main-sequence F7 stars. High inclination orbits can be 
excluded by the lack of eclipses from photometric moni- 
toring. Baliunas et al. (1997) exclude i ;> 83°. This is 
consistent with our experiment as we find no evidence for 
a companion at these high inclinations. 

We reiterate that the derivation of an upper limit for 
the geometric albedo requires the assumption of a value 
for the planetary radius (1.2 R p ) and a functional form 
for the phase variation (a Lambert law sphere). If the ac- 
tual values are significantly different than these, then the 
upper limit on the geometric albedo is modified as well. 
For example, assuming a smaller planetary radius would 
permit a larger albedo (see equation |l|) . 

Published predictions of the albedo of CEGPs vary by 
many orders of magnitude, and are highly sensitive to the 
presence of condensates in the planetary atmosphere. Bur- 
rows & Sharp (1999) consider cloud formation and de- 
pletion by rainout, and demonstrate that MgSiC>3 will 
be an abundant condensate at the effective temperature 
of t Boo b (~ 1500 K). Marley et al. (1999) calculate 
both cloud-free and silicate cloud atmospheres and pre- 
dict 0.35 <; p(480nm) < 0.55 for an EGP with a temper- 
ature of 1000 K, which is greater than our upper limit of 
p(480nm) = 0.3. They neglect the effects of stellar insola- 
tion on the model atmosphere. Seager & Sasselov (1998) 



explicitly include the stellar flux, solve the equation of ra- 
diative transfer through a model atmosphere of r Boo b, 
and predict p(480nm) ~ 0.0002. The low albedo is due 
in part to the absorption of photons by TiO in the atmo- 
sphere. However, it may be that the TiO forms and rains 
out, and thus is not an important factor. Including the 
presence of MgSiOs clouds, Seager & Sasselov predict a 
larger (but still very dark) albedo of p(480nm) ~ 0.003. 
The reflectivity of the MgSiOs grains at a given wave- 
length is highly dependent on the grain size relative to 
the wavelength of light. Burrows & Sharp (1999) also pre- 
dict that other condensates (such as Fe) may be present at 
these temperatures. If iron droplets are a significant con- 
densate, the resulting albedo would be very dark due to 
the high absorption at optical wavelengths. Given the cur- 
rent uncertainty in the models, there are many reasonable 
model planetary atmospheres which are consistent with 
our upper limit. 

We have achieved the current upper limit using only a 
limited spectral range, and data obtained when the planet 
was far from opposition. It is restricted by the photon 
noise of the data set, not by systematic errors. By ex- 
panding the spectral range and observing on several nights 
when the planet is near opposition, it would be possible to 
significantly reduce this upper limit. It may be advanta- 
geous to conduct this experiment at shorter wavelengths, 
since Seager & Sasselov (1998) predict a dramatic rise in 
the albedo shortwards of 420 nm. 
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Fig. 1. — The solid curves are the 90%, 99%, & 99.9% confidence levels on the upper limit for the relative reflected flux e as a function 
of orbital inclination i, if the reflected light is a copy of the stellar spectrum. The dashed curves are the same confidence levels under the 
assumption that the system is tidally locked and thus the planet reflects a non-rotationally-broadened copy of the stellar spectrum. Upper 
limits on the geometric albedo p under the assumption that R p = 1.2i?j up are shown on the right-hand axis, and the values for Jupiter and 
Uranus are included for comparison. The lack of transits excludes i > 83° , and i < 17° can be excluded under the assumption that the stellar 
rotation axis is co-aligned with that of orbital motion, as discussed in the text. 



